Unexpected Stoichiometries in Mg-O System under High Pressure 
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Using ab initio evolutionary simulations, we explore all the possible stoichiometries for Mg-0 
system at pressures up to 850 GPa. In addition to MgO, our calculations find that two extraordi- 
nary compounds MgC>2 and MgsC^ become thermodynamically stable at 116 GPa and 500 GPa, 
respectively Detailed chemical bonding analysis shows large charge transfer in all magnesium ox- 
ides. Mg02 contains peroxide ions [O-O] 2_ , while non-nuclear electron density maxima play the 
role of anions in the electride compound Mg3 02. The latter compound is calculated to have a much 
narrower band gap compared to MgO and Mg02- We discuss conditions at which Mg02 and Mg3 02 
might exist in planetary conditions. 



Magnesium oxide is one of the most abundant phases in 
planetary mantles, and understanding its high-pressure 
behavior is essential for constructing models of the 
Earth's and planetary interiors. For a long time, MgO 
was believed to be among the least polymorphic solids 
- only the NaCl-type structure has been observed in ex- 
periments at pressures up to 227 GPa pQ. Static theoret- 
ical calculations have proposed that the NaCl-type (Bl) 
MgO would transform into CsCl-type (B2) and the tran- 
sition pressure is approximately 490 GPa at K (474 GPa 
with the inclusion of zero-point vibration) [2]-[4]. Cal- 
culations also predicted that MgO remains non-metallic 
up to extremely high pressure (20.7 TPa) [3 j, making 
it to our knowledge the most difficult mineral to metal- 
ize. Thermodynamic equilibria in the Mg-0 system at 
0.1 MPa pressure have been summarized in a previous 
study [5 , concluding that only MgO is a stable composi- 
tion, though metastable compounds (Mg02, Mg04) can 
be prepared at a very high oxygen fugacities. 

Although MgO seems to be so simple and well under- 
stood, surprises may be in store for the researcher. For 
instance, it was shown recently that in the Li-H system in 
addition to the "normal" LiH, new counterintuitive com- 
pounds LiH2, LiHe and LiHg become stable at pressures 
above 100 GPa [6]. If similarly unusual stoichiometries 
become stable also in the Mg-0 system, this could have 
important planetological and chemical implications. In 
this paper, we explore this possibility and indeed find 
two new stoichiometries to be thermodynamically sta- 
ble at high pressures. These two new stable compounds 
(Mg02 and Mg302) exhibit interesting crystal structures 
with unusual chemical bonding and insulating and semi- 
conducting electronic structures, respectively. 

THERMODYNAMICALLY STABLE 
COMPOUNDS IN THE MG-O SYSTEM 

We have performed structure searches with up to 20 
atoms in the unit cell at pressures in the range 0-850 
GPa for the Mg-0 system. These searches yielded MgO 



as a stable oxide, but additionally two new compounds 
were predicted to be stable in the regions of high and low 
oxygen chemical potential, respectively. To confirm this 
and to obtain the most accurate results, we then focused 
search on two separate regions of chemical space; Mg- 
MgO and MgO-Mg, respectively. Since the structures 
in the two regions exhibit different properties, we report 
them separately. 



Exploring phase stability at high oxygen fugacities: 
MgO-O system and stable peroxide Mg02 

It is well known that monovalent (H-Cs) and divalent 
(Be-Ba and Zn-Hg) elements are able to form not only 
normal oxides, but also peroxides and even superoxides 
[7] (for instance, Ba02 has been well studied at both am- 
bient and high pressure 0[9]). Our structure prediction 
calculations did confirm that magnesium peroxide exists 
at ambient pressure, and has a cubic pyrite-type struc- 
ture, with Pa3 symmetry and 12 atoms in the unit cell - 
in good agreement with experimental results [10] . In this 
cubic phase, Mg is octahedrally coordinated by oxygen 
atoms (which form O2 dumbbells), see Fig. 2a. How- 
ever, Pa3 Mg02 (c-Mg02 from now on) is calculated to 
have a positive formation enthalpy from Mg and O2, and 
is therefore metastable. The calculation shows that on 
increasing pressure, c-Mg02 transforms into a tetragonal 
form with the space group 14/ mem. In the t-Mg02 phase 
(Fig. 2b), Mg is 8-coordinate. Here we see the same trend 
of change from 6-fold to 8-fold coordination as in the pre- 
dicted B1-B2 transition in MgO. but in Mg02 it happens 
at mere 53 GPa, compared to 490 GPa for MgO. Most 
remarkably, above 116 GPa the t-Mg02 structure has a 
negative enthalpy of formation from MgO and O2, indi- 
cating that t-Mg02 becomes thermodynamically stable. 
Furthermore, its stability is greatly enhanced by pressure 
and its enthalpy of formation becomes impressively nega- 
tive, -0.43 eV/atom, at 500 GPa ! We also examined the 
effect of temperature on its stability by performing quasi- 
harmonic free energy calculations. Inclusion of thermal 
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FIG. 1. (a) Convex hull for the MgO-0 system at high pres- 
sures; (b) the formation enthalpy as a function of pressure. 
For oxygen, the structures of the ("-phase [TT] and e-phase [12] 
were used. For MgO, Bl and B2 phases were considered, o 
Mg02 was considered at 50 GPa, while £-Mg02 was used at 
pressures higher than 50 GPa. 
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FIG. 2. Crystal structures of (a) c-Mg02 phase at 50 
GPa, space group Pa3, a=4.5236 A, Mg(0, 0, 0), 0(0.4074, 
0.4074, 0.4074); (b)£-Mg0 2 at 500 GPa, space group 74/ mem, 
a=3.3773 A, c=3.9849 A, Mg(0, 0, 0.75), 0(0.1260, 0.3740, 
0.5). The green polyhedra are drawn to show the coordination 
environment for the selected Mg atoms (6-fold in c-Mg02, and 
8-fold in t-Mg0 2 ). 



effects clearly does not bring any qualitative changes in 
the formation free energy (AG), and Mg02 remains sta- 
ble at high temperatures. 



Phase diagram of the Mg-MgO system: Mg3C>2 is an 
exotic stable compound 



the transition to fee and simple hexagonal phase at 456 
and 756 GPa, respectively. These results are in excellent 
agreement with previous studies [T3UT5] . Mg-rich oxides, 
such as Mg20 and Mg302 begin to show very competi- 
tive enthalpy of formation at pressures above 100 GPa. 
However, they are still not stable against decomposition 
into Mg and MgO, and their crystal structures could 
be thought as a combination of blocks of Mg and Bl- 
MgO. This situation qualitatively changes at 500 GPa, 
where we find that Mg302 becomes thermodynamically 
stable. This new stable phase has a very unusual tetrag- 
onal structure with the space group PA/mbm. 

This crystal structure can be viewed as a packing of O 
atoms and ID-columns of almost perfect body-centered 
Mg-cubes. As shown in Fig. 3, there are two types of 
Mg atoms in the unit cell, Mgl and Mg2 - Mg2 atoms 
form the cubes, filled by Mgl atoms. Within the cubic 
columns, one can notice empty (Mgl)2(Mg2)4 clusters 
with the shape of flattened octahedra, with Mg-Mg dis- 
tances ranging from 2.08 A (Mgl-Mg2) to 2.43 A (Mg2- 
Mg2). The coordination environments are quite differ- 
ent: each Mgl is bonded to two Mgl atoms and eight 
Mg2 atoms, while each Mg2 atoms is bonded to six O 
atoms (trigonal prismatic coordination). Oxygen atoms 
in £-Mg302 are coordinated by eight Mg2 atoms. 

What are the implications of these two compounds for 
planetary sciences? High pressures, required for their sta- 
bility, are within the range corresponding to deep plane- 
tary interiors. In planetary interiors, reducing conditions 
dominate, related to the excess of metallic iron. However, 
given the diversity of planetary bodies it is not impossi- 
ble to imagine that on some planets strongly oxidized 
environments can be present at depths corresponding to 
the pressure of 116 GPa (in the Earth this corresponds 
to depths below 2600 km). At the more usual reduc- 
ing conditions of planetary interiors, Mg 3 02 could ex- 
ist at pressures above 500 GPa in deep interiors of gi- 
ant planets. There, it can coexist in equilibrium with 
Fe (but probably not with FeO, according to our DFT 
and DFT+U calculations of the reaction Fe + 3MgO — » 
FeO + Mg302). According to our calculations (Fig. 4), 
Mg302 can only be stable at temperatures below 1800 
K, which is too cold for deep interiors of giant planets; 
however, impurities and entropy effects stemming from 
defects and disorder might extend its stability field into 
planetary temperatures. Exotic compounds Mg02 and 
Mg 3 02, in addition to their general chemical interest, 
might be important planet-forming minerals in deep in- 
teriors of some planets. 




For the Mg-rich part of Mg-0 phase diagram, USPEX 
shows completely unexpected results. First of all, ele- 
mental Mg is predicted to undergo several phase tran- 
sitions induced by pressure: hep - bec - fee - sh. At 
ambient condition, Mg adopts the hep structure, while 
bcc-Mg is stable from 50 GPa to 456 GPa, followed by 



ELECTRONIC STRUCTURE 

Phonon calculations for Mg302 and Mg02 at pres- 
sures of their stability show that no imaginary phonon 
frequencies exist throughout Brillouin zone, suggesting 
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FIG. 3. Crystal structures of £-Mg302 at 500 GPa, space 
group PA/mbm, a=4.5079 A, c=2.3668 A, Mgl(0.3494, 
0.1506, 0.5); Mg2(0, 0, 0); 0(0.8468, 0.6532, 0). This rep- 
resentation highlights the Mg octahedra. 
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FIG. 4. (a) Convex hull for the Mg-MgO system at high 
pressures. For Mg, the bcc, fee and sh structures were used 
at their stability range. For MgO, Bl and B2 phases were 
considered. £-Mg 3 02 structure was used in the whole investi- 
gated pressure range, (b) Predicted P — T phase diagram for 
the Mg-MgO system. 



that these structures are dynamically stable. Together 
with our calculated thermodynamic functions, this sug- 
gests thermodynamic stability of these compounds. How 
are these exotic magnesium oxides stabilized at high pres- 
sure? To answer this question, we analyzed the electronic 
structure and chemical bonding for these compounds. 

Electron Localization Function (ELF) gives informa- 
tion about the bonding character and valence electron 
configuration of atoms in a compound [T6] . The ELF 
pictures show asymmetry of the Mg-0 bonds, indicating 
significantly ionic character of bonding. As shown in Fig. 
[5^i, valence electrons in t-Mg02 are mainly concentrated 



around O atoms. Charge transfer was also investigated 
on the basis of the electron density using Bader analy- 
sis [HI [18]. In t-Mg0 2 , the net charge on Mg is +1.747 
e, indicating the nearly complete transfer of Mg 3s elec- 
trons to O atoms (just like Mg in MgO: Bader charges 
are +1.737 e at GPa and +1.675 e at 600 GPa, respec- 
tively). Each O has almost 7 valence electrons (6.873 
e), thus with the formation of singly bonded 0-0 dumb- 
bell the octet rule is fulfilled. Each 0-0 dumbbell can 
be viewed as a peroxide-ion [O-O] 2_ with a closed-shell 
electronic configuration. 



The ELF distribution in £-Mg 3 02 (Fig. |5|d) also con- 
firms strong charge transfer from Mg to O. However, we 
surprisingly found a very strong interstitial ELF maxi- 
mum (ELF=0.97) located in the center of Mg-octahedron 
(Fig. [5J3). To obtain more insight, we performed Bader 
analysis. The resulting charges are +1.592 e for Mgl, 
+1.687 e for Mg2, -1.817 e for O, and -1.311 e for the 
interstitial electron density maximum. Such a strong 
interstitial electronic accumulation requires an explana- 
tion. Its electronic structure (Fig. |5^) exhibits intrigu- 
ingly high occupancy of Mg-p, d-orbitals with overlap- 
ping energy ranges, which implies a strong Mg 3p-3d 
hybridization. At high pressure, strong interstitial elec- 
tron localization was found in alkali and alkaline-earth 
elements, for instance, sodium becomes insulator due to 
strong core-core orbital overlap [19]. As a measure of 
size of the core region we use the Mg 2+ ionic radius (0.72 
A [20]), while the size of the valence electronic cloud is 
represented by the 3s orbital radius (1.28 A [21 J. In 
Mg 3 2 , Mg-Mg contacts at 500 GPa (2.08 A for Mgl- 
Mg2, 2.37 A for Mgl-Mgl and 2.43 A for Mg2-Mg2) are 
only slightly shorter than the sum of valence orbital radii. 
Thus, the main reason for strong interstitial electronic 
localization is the formation of strong multicenter cova- 
lent bonds between Mg atoms; the core- valence expulsion 
(which begins at distances slightly longer than the sum 
of core and valence radii and increases as the distance 
decreases) could also be an important factor for valence 
electron localization. Strong Mg-Mg covalent bonding is 
not normally observed; the valence shell of the Mg atom 
only has a filled 3s 2 orbital, a configuration unfavorable 
for strong bonding. Under pressure, the electronic struc- 
ture of the Mg atom changes (p- and d-levels become sig- 
nificantly populated), and strong covalent bonding can 
appear as a result of the p-d hybridization. There is an- 
other way to describe chemical bonding in this unusual 
compound. We remind that Mg 3 02 is anion- deficient 
compared with MgO; the extra localized electrons in Mg 
octahedron interstitial play the role of anions, screening 
Mg atoms from each other. These two descriptions are 
complementary. 
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FIG. 5. ELF isosurface of a) t-Mg0 2 structure; b) t-Mg 3 2 structure at 500 GPa and ELF = 0.83. c) Charge density 
distribution of £-Mg3 02 along c axis. Band structure and partial densities of states for d) t-tetragonal MgC^ and e) £-Mg3 02- 
f) Band gaps as a function of pressures for various Mg-0 compounds. 



CONCLUSIONS 



METHODS: CRYSTAL STRUCTURE SEARCH 



In summary, we performed a systematic search for pos- 
sible stoichiometrics in the Mg-0 system at pressures 
up to 850 GPa. Other than the well-known compound, 
MgO, we found that two more stoichiometrics (Mg02 and 
Mg302) become thermodynamically stable at pressures 
above 116 GPa and 500 GPa, respectively. Our analy- 
sis reveals that bonding in both t-Mg02 and t-Mg 3 02 
exhibits significantly ionic character. Mg02 is stabilized 
by the formation of the peroxide [O-O] 2_ anion, while the 
strong electron localization in Mg octahedron plays the 
role of an additional anion and makes the anion-deficient 
compound Mg 3 02 stable. Compared with MgO, which is 
a wide gap insulator (DFT band gap 7.73 eV at 500 GPa 
- see Fig. |5]F), Mg302 is predicted to have a much nar- 
rower band gap (it is a semiconductor with a DFT gap of 
2.20 eV), while Mg0 2 is an insulator with the DFT band 
gap of 5.35 eV; note that DFT gaps are always lower 
bound to the true gaps. These two compounds might ex- 
ist in interiors of other planets, and are thus potentially 
important for our fundamental understanding of the uni- 
verse - in addition to presenting striking new chemical 
phenomena. 



Searches for the stable compounds and structures 
were performed using evolutionary algorithm USPEX, 
as implemented in the USPEX code[22j [23]. This was 
done in conjunction with ab initio structure relaxations 
based on density functional theory (DFT) within the 
Perdew-Burke-Ernzerhof (PBE) generalized gradient ap- 
proximation (GGA)[24], as implemented in the VASP 
code [25]. For structural relaxation, we used the all- 
electron projector- augmented wave (PAW) method [26] 
and plane wave basis set with the 600 eV kinetic energy 
cutoff; Brillouin zone was sampled by Monkhorst-Pack 
meshes with the resolution 2tt x 0.06 A -1 . Such cal- 
culations provide an excellent description of the known 
structures (Mg, O2, MgO) and their energetics. The 
most significant feature of USPEX we used in this work 
is the capability of searching for a specific area of the 
composition space - as opposed to the more usual struc- 
ture predictions at fixed chemical composition. The de- 
sired composition space is described via building blocks 
(for example, search for all compositions in a form of 
[xA^Oa+^/MgO] or [xMg+yO]). During the initializa- 
tion USPEX samples the whole range of compositions 
of interest randomly and sparsely. Chemistry-preserving 
constraints in the variation operators are lifted and re- 
placed by block-correction scheme which ensures that a 
child structure is within the desired area of composition 
space. A new "chemical transmutation" operator is in- 
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troduced [27J [28]. Stable compositions are determined 
using the convex hull construction: a compound is ther- 
modynamically stable if the enthalpy of its decomposi- 
tion into any other compounds is positive. To ensure 
that the obtained structures are dynamically stable, we 
calculated phonon frequencies across the Brillouin zone 
using the finite-displacement approach as implemented 
in the Phonopy code [29]. We also explored the effects 
of temperature using the quasiharmonic approximation; 
for each structure, phonons were calculated at 10 - 15 
different volumes. 

Calculations were performed at the supercomputer of 
Center for Functional Nanomaterials, Brookhaven Na- 
tional Laboratory, which is supported by the U.S. De- 
partment of Energy, Office of Basic Energy Sciences, un- 
der contract No. DE-AC02-98CH10086, at the Joint 
Supercomputer Center of the Russian Academy of Sci- 
ences and on Skif supercomputer (Moscow State Uni- 
versity). This work is funded by DARPA (grants 
No. N660011014037, W31P4Q1210008), National Sci- 
ence Foundation (No. EAR-1114313). 
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